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Abstract 

We propose a new analytical method to solve for nonexactly soluble 
Schrodinger equation via expansions through some existing quantum 
numbers. Successfully, it is applied to the rational nonpolynomial 
oscillator potential. Moreover, a conclusion reached by Scherrer et al. 
[2], via matrix continued fractions method, that the shifted large N 
expansion method leads to dubious accuracies is investigated. The 
cutoff - Coulomb and Coulomb plus logarithmic potentials are also 
investigated. 



1 



1 Introduction 



The importance of rational nonpolynomial oscillator (NPO) potential, 

V ^ = a °v 2 + rfw ' 6>0 ' (1) 

arises in nonlinear Lagrangian field theory, laser theory, and elementary par- 
ticle physics [1-4]. Only a class of exact analytical solutions for a certain 
parameter dependence a = a(b) obtains [3,5-7]. Hence, it has been a subject 
of several investigations [1-15] ( exhaustive lists of these could be found in 
Ref.s[2,15]). 

On the other hand, the cutoff - Coulomb potential 

V(q) = ]- ,c>0, (2) 

q + c 

with the truncation parameter c, avoids the singularity at q = in the 
Coulomb potential ( the crux of divergence difficulties in quantum field the- 
ory [16-18]). It has been suggested [18] that if gravitational interactions of 
elementary particles are taken into acount, there would be a gravitational 
cutoff of Coulomb interaction. Equation (2) represents a nonrelativistic ex- 
pression of this idea. It also serves as an approximation to the potential due 
smeared charge rather than a point charge. 

The Schrodinger equation with such interactions, (1) and (2), belongs to 
the class of quantum mechanical systems which are nonexactly soluble in 
general. However, the solutions of exactly solvable potentials ( an interesting 
field of mathematical physics in itself [19,20]) can be used in perturbation 
and pseudoperturbation theories, or they can be combined with numerical 
calculations. Nevertheless, in the simplest case, analytical calculations can 
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aid numerical studies in areas where numerical techniques might not be safely 
controlled. For example, when bound - state wave functions with arbitrary 
nodal zeros are required for certain singular potentials (a next level of com- 
plexity), analytical solutions can supply a basis for numerical calculations. 
Moreover, in many problems the Hamiltonian does not contain any physi- 
cal parameter suitable for a perturbation expansion treatment. More often, 
the Hamiltonian contains physical parameters, but, typically, zeroth - order 
solutions for special values of these are not tractable or good starting approx- 
imations. Alternatively, one would resort to apparently artificial conversions 
to perturbation problems which have been shown to make progress [21-35]. 

The three - dimensional (3D) spherically symmetric NPO potential has 
been investigated [8] by means of the shifted large - N expansion technique 
(SLNT). However, Scherrer and co-workers [2] have employed a matrix con- 
tinued fractions (MCF) method and concluded that SLNT [8] leads to dubi- 
ously accurate results in the critical range a 1 / 2 jb ps 0.1 — 30. Handy et al. 
[14] have used the eigenvalue moment method (EMM) to obtain upper and 
lower energy bounds, for I ^ 0, and compared their results with those of Roy 
et al. [12] by the suppersymmetric quantum mechanics (SSQM). Singh et al 
[17] have used a numerical integration method (NIM) to find the eigenvalues 
for the 3D cutoff - Coulomb potential. Hall and Saad [35] have used a smooth 
transformation (STM) and a numerical integration methods to obtain bound 
- states for the Coulomb plus a logarithmic perturbation term. Hence, apart 
from those of SLNT, results from convincingly powerful methods exist for 
comparison purposes. 

Recently, we have introduced a pseudoperturbative shifted - 1 ( I is the an- 
gular momentum quantum number) expansion technique ( PSLET) to solve 
for nodeless states of Schrodinger equation. It simply consists of using 1/7 as 
a pseudoperturbation parameter, where I = I— (5 and (3 is a suitable shift. The 
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shift (3 is vital for it removes the poles that would emerge, at lowest orbital 
states with 1=0, in our proposed expansions below. Our analytical, or often 
semianalytical, methodical proposal PSLET has been successfully applied to 
quasi - relativistic harmonic oscillator [34], spiked harmonic oscillator [32], 
and to anharmonic oscillators [35] potentials. 

Encouraged by its satisfactory performance in handling nodeless states, 
we generalize PSLET recipe ( in section 2) for states with arbitrary number 
of nodal zeros, k > 0. Moreover, in the underlying "radical" time - indepen- 
dent radial Schrodinger equation, in h — m — 1 units, 



1<P_ 1(1 + 
~2~cLf + 2q 2 



1) 



+ V{q) 



(3) 



the isomorphism between orbital angular momentum I and dimensionality 
D invites interdimensional degeneracies to obtain [36-39] (for more details 
the reader may refer to ref.s [36,38]). Hence, the symmetry of an attendant 
problem obviously manifests the admissibility of the quantum number /: In 
one-dimension (ID), / specifies parity, (— with the permissible values 
-1 and/or ( even and/or odd parity, respectively) where q = x G (— oo, oo). 
For two-dimensional (2D) cylindrically symmetric Schrodinger equation one 
sets I = \m\ — 1/2, where m is the magnetic quantum number and q = (x 2 + 
y 2 ) 1 / 2 g (0, oo). Finally, for three-dimensional (3D) spherically symmetric 
Schrodinger equation, / denotes the angular momentum quantum number 
with q = (x 2 + y 2 + z 2 ) 1 ' 2 G (0, oo). 

In section 3, we investigate PSLET recipe and consider, for the sake of 
diversity, the potentials; (i) V(q) = A 2 q 2 /2 : the harmonic osillator ( the limit 
of (1) when b — > oo and/or a — > 0, with a Q = A 2 /2), (ii) V(q) = —l/q, 
the Coulomb ( the limit of (2) when c — ► 0, (iii) the NPO (2), (iv) the 
cutoff - Coulomb (2), and (v) the Coulomb perturbed by a logarithmic term, 
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which has no experimental evidence, to the best of our knowledge, thus our 
calculations are only of academic interest. The last section is reserved for 
summary and remarks. 



2 The generalization of PSLET 

The construction of our PSLET starts with shifting the angular momentum 
in (3) to obtain 

[ Id 2 I 2 + (2/3 + 1)1 + f3(/3 + 1) P . .1 _ .. _ _ . . fA . 

where Q is a constant that scales the potential V(q) at large - I limit ( the 
pseudoclassical limit [36]) and is set, for any specific choice of / and k, equal 
to I 2 at the end of the calculations. Next, we shift the origin of the coordinate 
system through x = l l ^ 2 (q — q )/q , where q Q is currently an arbitrary point 
to be determined below. Expansions about this point, x = (i.e. q = q ), 
yield 

^ = E(-ir ^ r n! \ (5) 

« n=0 % 



v(x(q)) = ± (*m) is^L r~». (6) 

Obviously, the expansions in (5) and (6) localize the problem at an arbitrary 
point q Q and the derivatives, in effect, contain information not only at q Q but 
also at any point on g-axis, in accordance with Taylor's theorem. It is then 
convenient to expand E^i as 
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M i-n 



E k ,i - J2 E Z 1 



(7) 



n=-2 



Equation (4) thus becomes 



(8) 



with 



jV{x{q)) 



1 V(go) 

2g 2 Q 



+ (2/3 + 1) f;(_i)»ll!±i x »rn/2 + f; B n x n r^ 2 



n=l 

oo 



n=3 



+ / g(/5 + i)^(-if^^x"r( n+2 )/ 2 , 



(9) 



n=0 



n+2 
o 



dq™ J n\Q 



(10) 



Equation (8), along with (9) and (10), is evidently the one - dimensional 
Schrodinger equation for a perturbed harmonic oscillator 



X k (x) = \ k X k (x), 



(11) 



where w 2 = 2B 2 , 



e Q = I 



1 QqV(Qq) 

2 Q 



2/3+1 13(13 + 1) 
2 21 



(12) 
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and P(x) represents the remaining terms in eq.(9) as infinite power 
perturbations to the harmonic oscillator. One would then imply that 



Xh — I 



1 

+ 7 



1 ffl 

2 Q 

W+l) X (Q) 

o r- \ 



2/3 + 1 n 1. 

oo 

+ E 4 n_1) ^ n , 

n=2 



and 



n=-2 



Hence, equations (13) and (14) yield 



F (-2)_ 1 , ^(9o) 

" 2<g + Q 



E 



(-i) 



2/3 + 1 1. 



E {0) - — 



W+l) + A (o) 



4? = A^o 2 ; n > 1. 



Where g c is chosen to minimize E^^, i. e 



dE, 



(-2) 



dq Q 



= and 



d 2 E, 



(-2) 
k,l 



> 0. 



(19) 



Hereby, V(q) is assumed to be well behaved so that e\, ^ has a minimum q Q 
and there are well - defined bound - states. Equation (19) in turn gives, with 

i=VQ, 

(20) 



l-P = MV(q ). 



Consequently, the second term in Eq.(9) vanishes and the first term adds a 
constant to the energy eigenvalues. It should be noted that the energy term 
P-E^ 2 - 1 corresponds roughly to the energy of a classical particle with angular 
momentum L z =l executing circular motion of radius q in the potential V(q ). 
It thus identifies the zeroth - order approximation, to all eigenvalues, as a 
classical approximation and the higher - order corrections as quantum fluc- 
tuations around the minimum q , organized in inverse powers of I. The next 
correction to the energy series, lE^, consists of a constant term and the 
exact eigenvalues of the harmonic oscillator w 2 x 2 /2.The shifting parameter 
(3 is determined by choosing TE^^O. This choice is physically motivated. 
In addition to its vital role in removing the singularity at I — 0, it also re- 
quires the agreements between PSLET eigenvalues and eigenfunctions with 
the exact well known ones for the harmonic oscillator and Coulomb poten- 
tials. Hence 



1 n l ^ ' 

2 + ( ' + 2> 



w = 



3 | g° V "fa>) (21) 



where primes of V(q ) denote derivatives with respect to q Q . Then equation 
(9) reduces to 
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Q fV(x(q)) = q 2 J 



1 V(qo 

2g 2 Q 



n=0 



(22) 



where 



v^(x) = B 2 x 2 + ^±A, 



(23) 



(24) 



t; (n) (x) = £ n+2 x n+2 + (-1)™ (2/? + 1) 



n+2 



(ra + 1) 



+ (_l)n W+l) ( n _ 1) > 



n > 2. 



(25) 



Equation (8) thus becomes 



1 J2 oo 



oo 

Qo^k,i 1 



.71=1 



(26) 



Up to this point, one would conclude that the above procedure is nothing 
but an imitation of the eminent shifted large-N expansion (SLNT) [24-29,35]. 
However, because of the limited capability of SLNT in handling large-order 
corrections via the standard Rayleigh-Schrodinger perturbation theory, only 
low-order corrections have been reported, sacrificing in effect its preciseness. 
Therefore, one should seek for an alternative and proceed by setting the wave 
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functions with any number of nodes as 



tffc,j(a;(g)) = F k)i (x) exp(U k ,i(x)) . 



(27) 



In turn, equation (26) readily transforms into the following Riccati equation: 



1 oo 

■« (^» + ^(^iiW) + E « (B) (*)i- 



n/2 



n=0 



E^r^ 



n=l 



(28) 



where the primes denote derivatives with respect to x. It is evident that this 
equation admits solution of the form 



u'^x) = E ut\x) r n/2 + £ G[ n \x) r«+w (29) 



?1=0 



n=0 



00 fe— 1 



n=0 p=0 



(30) 



where 



n+1 



P) 2m-i . n — n 

/ y J - y m,n,fe x 1 - i - y 0,n,fe u ; 

m=0 



(31) 



71+1 



^aT ( a; ) — E 

m=0 



(32) 



Substituting equations (29) - (32) into equation (28) implies 
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("[/( m )f/( Tt - m )/-™/2 _|_ gHg("-"»)p(n+2)/2 
/ 4 \ k k k k 

n=0 m=0 



+ 2ut ] Gt~ m) r {n+1)/2 ) + E v {n) r n / 2 - E <z 2 4?f ^ 



—n 



n=0 n=l 



00 1 

^(*) E (Ui n) r^ + Gpn-W) -^» = (33) 

.71=0 J Z 



The above procedure obviously reduces to the one described by Mustafa and 
Odeh [33,34,36,37], for k — 0. Moreover, the solution of equation (33) follows 
from the uniqueness of power series representation. Therefore, for a given k 
we equate the coefficients of the same powers of I and x, respectively. For 
example, when k — 1 one obtains 



£1,0,1 = -w, 




= — W X, 



(34) 



B 3 




0,0,1 



(35) 



w 



w 



1(2^,0,1 + 2/3 + 1), 



(36) 




(37) 
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A,2,i = 1(1 £2,2,1 + Co,o,i Ci,o,i - |(2^ + 1)) , (38) 




3 Dj 
2 



1,2,1 



CI 



l f) , (39) 



E 



■o,i ^1,0,1 - 



etc. Here, we reported the nonzero coefficients only. One can then calcu- 
late the energy eigenvalues and eigenfunctions from the knowledge of C m , n ,k-i 
D m n k, an d a£l in a hierarchical manner. Nevertheless, the procedure just 
described is suitable for a software package such as MAPLE to determine 
the energy eigenvalue and eigenfunction corrections up to any order of the 
pseudoperturbation series, (7) and (29)-(30). 

Although the energy series, equation (7), could appear divergent, or, at 
best, asymptotic for small I, one can still calculate the eigenenergies to a very 
good accuracy by forming the sophisticated [N,M] Pade approximation [21] 




(40) 



by 



and with the P\ 



(1//) Pade approximant it becomes 



E k ,[AA]=PEif+P 4 \l/l). 



(41) 



Our recipe is therefore well prescribed. 
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3 Some applications 



We begin with the harmonic oscillator potential V(q) = A 2 q 2 /2 and find q Q 
from (20), along with (21). Once q Q is determined, Eq.(15) reads 



E ( ~ 2) P = Al ; I = 2k + 1 + 3/2, (42) 

Ei?l = Ef) = E^ = Ef]T 2 = ... = ^ (43) 

Hence, the corresponding energies are 

Eo,i = A (l + |) , (44) 

the well known exact results. For k — 



tt t \ l ( y \ y3 y \ y5 y \ \ 
u ^ x) = -2{ y -Y + j-^ + j--6 + ---) 

,y( y 2 . y 3 y\ y 5 y 6 . \ 

+/ [y-Y + Y-^ + j-j + ---) 

-l\-Ty ; y = xi- 1 ' 2 . (45) 



Obviously the terms between brackets in equation (44) are the infinite geo- 
metric series expansions for ln(l + y). Equation (44) can be recast as 

U 0>l (x) = ln(l + y)- 1 ' 2 + ln(l + y) T -ly - f-, (46) 
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which in turn implies the exact wave functions 

*o,i(?) = Nojf-We-**' 2 ; I = Aq 2 , (47) 

where N j are the normalization constants. 

Next we consider the Coulomb potential V(q) — —1/q. In this case, 
Eq.(15) reads 



ElT^ = -^2 S l = k + l + l = ql (18) 



The reminder corrections are identically zeros, Eq.(43). Hence, one obtains 
the eigenvalues 

Ek > 1 = 2(k+i+ir m 

the well known exact energies for the Coulomb potential. For k = 

U ,i(x) = -ly + ln(l + y) T , (50) 

which in turn implies the exact wave functions 

*oM= N <>We- iq/qo - ( 51 ) 
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Now let us consider the NPO potential (1), for which Eq.(21) implies 



/3 = ~[l + (2k + l)w] ; 



w 



2 



(1 + 3bg 2 Q + 3b 2 qt + 6 3 gg + a) 
\(l + 2bq* + Pq* + a)(l + bqiy 



(52) 



In turn Eq.(20) reads 




1 + bql 



^l + a + bql + b 2 qt. 



(53) 



Equation (43) is explicit in q Q and evidently a closed form solution for q 
is hard to find, though almost impossible. However, numerical solutions are 
feasible. Once q Q is determined the coefficients C m;n and D m n are obtained in 
a sequel manner. Consequently, the eigenvalues, Eq.(40), and eigenfunctions, 
Eqs.(29)-(32), are calculated in the same batch for each value of a, b, and I. 

In order to make remediable analysis of our results we have calculated the 
first ten terms of the energy series. The effect of each term has been taken 
into account. We have also computed the Pade approximants E[N,M] for 
N = 2, 3,4 and M = 2,3,4. Therefore, the stability of the sequence of the 
Pade approximants was in point. 

In tables 1 and 2 we list PSLET results E P , Eq.(40), along with the [4,4] 
Pade approximants, Eqs.(41). The results of Roy et al.[12], via SSQM, and 
Handy et al. [14] , via EMM, are also displayed for comparison. In tables 3-5 
we compare our results with those of Scherrer et al.[2], via MCF, and Varshni 
[8], via SLNT. 

In tables 1 and 2 we have observed that, if the last two digits are neglected, 
the energy series Eq.(40) stabilizes at E 5 up to E w . Where E 5 and E w denote 
that the energy is computed by the first five and first ten terms of the energy 
series, respectively. Yet for Z = 5 — 20 it stabilizes at E 3 up to E w . While the 
Pade approximants for / = 5 — 20 had no effect on the energy series their effect 
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is not dramatic for / = 1 — 3. Also it should be reported that the sequence 
of Pade approximants stabilizes at E[3, 3] for 1 — 1 — 3 and at E[2, 2] for 
/ = 5 — 20. Therefore, one could confidently conclude that the results from 
the Pade approximants are exact provided that some uncertainty lies in the 
last two digits for 1 — 1 — 3. Nevertheless, our results do not contradict with 
the upper and lower bounds computed from EMM [14]. To a satisfactory 
extent they also agree with those from SSQM [12]. 

In tables 3 - 5 we list our results along with those of Scherrer [2], from 
MCF, and Varshni [8], from SLNT. Our results compare better with those 
of Scherrer than the results of SLNT. However, collecting only the first four 
terms of the energy series Eq.(40) we found that PSLET results are in exact 
accord with SLNT. Moreover, the trends of convergence of the energy series 
and the sequence of Pade approximants are similar to those for / = 1 — 3 in 
tables 1 and 2. 

In table 6 we display our results along with those from STM and direct 
numerical integration (DNI) reported in [35] for the Coulomb plus logarithmic 
potential. In table 7 we list our predictions for the cutoff - Coulomb and 
compare them with those of Singh et al. [17]. Finally, we report the k=0 
eigenvalues for the NPO and cutoff - Coulomb potentials in table 8 and 9, 
respectively. 

4 Summary and Remarks 

In this work we have introduced a pseudoperturbative shifted - / expansion 
technique (PSLET) to deal with the calculation of the energy eigenvalues and 
eigenfunctions of Schrodinger equation in one batch. We have shown that it 
is an easy task to implement PSLET without having to worry about ranges of 
couplings and forms of perturbations in the potential involved. In contrast to 
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the textbook Rayleigh - Schrodinger perturbation theory, an easy feasibility 
of computation of the eigenvalues and eigenfuctions has been demonstrated, 
and satisfactory accuracies have been obtained. Perhaps it should be noted 
that for each entry in tables 1-5 one can construct the wavefunction from the 
knowledge of C m>n and D m ^ n . However, such a study lies beyond the scope 
of our methodical proposal. 

The conclusion reached by Scherrer [2], via MCF method for the NPO 
potential Eq.(l), that SLNT [8] leads to dubiously accurate results in the 
critical range a l l 2 /b 0.1 — 30 has been confirmed in the present work. 
However, the natural extensions of SLNT or its variants [8,16-21], repre- 
sented by PSLET or the modified SLNT [22,23], show that one could get 
convincingly reliable results. Moreover, the dubious accuracies of SLNT [8] 
should be attributed mainly to the limited capability of SLNT to calculate 
the energy series beyond the fourth - order term. Our results for the cutoff - 
Coulomb and Coulomb plus logarithmic potentials are also very satisfactory. 

A final remark is in point. The attendant method could be applied to 
systems at lower dimensions. Here is the recipe. Rewrite the centrifugal term 
in Eq.(2) as A(A + l)/2g 2 . Shift A through A = A — (3 and expand in inverse 
powers of A following exactly the same procedure described in section 2. In 
this case, A = I in three dimensions ( with q > 0), A = m — 1/2 in two 
dimensions where m is the magnetic quantum number ( with q > 0), and 
A = — 1 and/or in one dimension ( with — oo < q < oo) [23-28]. 
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Table 1: Bound - state energies of the NPO potential (1). Where E P 
represents PSLET results, Eq.(51), E ss from SSQM [12], E E from EMM [14]. 
The Pade approximant E[3, 4] and the upper bound of EMM are obtained, 
respectively, by replacing the last j digits of E[3, 3] and of the lower bound 
of EMM with the j digits in parentheses. 



/ 


a 


b 


Ep 


£[3,3] & (^[3,4]) 


Ess 


E E 


1 


0.1 


0.1 


5.1863731 


5.1863730 (29) 


5.186338 


5.1863730 (30) 






0.5 


5.100883 


5.100858 (61) 


5.100976 


5.100842 (65) 






1 


5.065556 


5.065563 (71) 


5.065610 


5.06428 (609) 




0.5 


0.1 


5.8935959 


5.8935951 (27) 


5.893494 


5.8935952 (52) 




1 


0.1 


6.7042393 


6.704239 (46) 


6.704090 


6.7042389 (89) 






1 


r /;ri no 


o.Dollo (4UJ 


r o rroi 1 O 

o.oozliz 


o.ooOo [Zl ) 




10 


0.1 


15.8137089 


15.81370943 (48) 


15.813628 


15.81370943 (43) 




100 


0.1 


49.38979430 


49.389794297 (97) 


49.389615 


49.38979427 (34) 






10 


14.371 


14.3622 (64) 


14.363739 


1.7 (14.609) 






100 


5.9934347 


5.993450 (59) 


5.993565 


— (6.389) 


2 


0.1 


0.1 


7.24396166 


7.243961847 (50) 


7.243927 


7.243961840 (40) 






0.5 


7.118983 


7.1189816 (10) 


7.119005 


7.11890 (901) 






1 


7.0737228 


7.0737258 (70) 


7.073713 


7.0730 (44) 




0.5 


0.1 


8.17787177 


8.17787168 (66) 


8.177754 


8.17787169 (69) 




1 


0.1 


9.2619150 


9.26191476 (42) 


9.261812 


9.26191478 (78) 






1 


7.73479 


7.734821 (33) 


7.734778 


7.734 (36) 




10 


0.1 


21.83609251 


21.83609251 (58) 


21.836043 


21.83609247 (54) 




100 


0.1 


68.802061155 


68.802061155 (55) 


68.801562 


68.8020606 (15) 






10 


16.61083 


16.61081 (54) 


16.611028 


16.5997 (6540) 






100 


7.9960247 


7.99602475 (91) 


7.996048 


7.9947 (8.0378) 
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Table 2: Same as table 1. 



1 


a 


b 


Ep 


£[3,3] & (£[3,4]) 


Ee 


3 


0.1 


0.1 


9.2943590 


9.294359116 (16) 


9.2943591109 (09) 






0.5 


9.1318120 


9.131807 (12) 


9.131799 (838) 






1 


9.0789113 


9.0789116 (13) 


9.0787 (92) 




0.5 


0.1 


10.4292039 


10.429204129 (33) 


10.42920412 (12) 




1 


0.1 


11.7606210 


11.760620955 (50) 


11.76062096 (96) 






1 


9.787665 


9.787668 (73) 


9.7875 (81) 




10 


0.1 


27.68830286 


27.6883029 (31) 


27.6883028 (28) 




100 


0.1 


88.01806590 


88.018065906 (07) 


88.0180658 (60) 






10 


18.719989 


18.720001 (21) 


18.7186 (307) 






100 


9.99715366 


9.997153642 (55) 


9.9969 (10.0113) 


5 


200 


0.1 


179.48311321 


179.48311321 (21) 


179.483107 (16) 




500 


0.1 


286.13076490 


286.130764905 (05) 


286.13073 (82) 




1000 


0.2 


401.608033488 


401.6080334886 (89) 


401.6078 (81) 




10000 


0.4 


1280.6255249 


1280.62552494 (95) 


1280.6254 (56) 




10000 


0.5 


1275.7839677 


1275.7839677 (77) 


1275.7838 (42) 


10 


200 


0.1 


311.86088089 


311.86088089 (89) 


311.8601371 (16266) 




1000 


0.1 


713.36153440 


713.36153440 (40) 


713.36081 (321) 




1000 


0.2 


699.10562257 


699.10562257 (57) 


699.10424 (909) 




10000 


0.4 


2242.79417589 


2242.79417589 (89) 


2242.7891995 (3867) 




10000 


0.5 


2228.5184345 


2228.5184345 (45) 


2228.513255 (30746) 


20 


500 


0.1 


914.36631099 


914.36631099 (09) 


914.36540 (851) 




1000 


0.1 


1312.251674809 


1312.251674809 (09) 


1312.25006 (333) 
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Table 3: Bound - state energies of the NPO potential (1) for a = 10. Ep, 
£[3,3], and £[3,4] are the same as in table 1. E M from MCF [2] and E SL 
from SLNT [8]. 



/ 

I 


h 
u 


E p 


J-/\0]0\ 06 1X^10,4:11 


Em 







-I 

1 


7.41837 


7.417532 (12) 


f-j 1 i Tr A/" 1 

7.417506 


7.4056 




10 


o o o r 

3.885 


3.864 (75) 


o o Tn no'7 

3.879037 


3.8732 




100 


3.0982o 


3.097906 (26) 


3.089317 


3.0984 




1000 


3.0099813 


o aaaa*7aa /oai\ 

3.0099799 (801) 


3.009981 


3.0100 


1 


-I 

1 


11.0714 


11.07326 (27) 


11.073300 


1 1 ATI /I 

11.0714 




1 a 

10 


r A /I non 

5.94089 


5.94057 (23) 


r (\ A AA^n 

5.940860 


r A A AO 

5.9408 




1 aa 

100 


r AAAO /I /I 

5.099344 


5.0993455 (64) 


r AAAQ /I /I 

5.099344 


r AAA A 

5.0994 




i nnn 

1UUU 


*J . \J \j ZJ ZJ ZJ o o o L ± 


^ nn999^^4fi f41 


O .\j\jZJZJZJO 


^ m nn 


2 


1 


14.0861 


14.08540 (31) 


14.085383 


14.0900 




10 


7.9622256 


7.962228 (57) 


7.962230 


7.9622 




100 


7.09960269 


7.099602614 (30) 


7.099603 


7.0996 




1000 


7.009996003 


7.00999600266 (38) 


7.009996 


7.0100 


3 


1 


16.719284 


16.71910 (33) 


16.719332 


16.7200 




10 


9.9724544 


9.9724556 (62) 


9.972455 


9.9724 




100 


9.09971542 


9.09971542 (42) 


9.099715 


9.0998 




1000 


9.00999714406 


9.00999714399 (92) 


9.009997 


9.0100 


4 


1 


19.137789 


19.137816 (24) 


19.137821 


19.1376 




10 


11.97836435 


11.97836459 (62) 


11.978365 


11.9784 




100 


11.099778408 


11.099778407 (07) 


11.099778 


11.0998 




1000 


11.00999777842 


11.009997778412 (64) 


11.009998 


11.0100 
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Table 4: Same as table 3 for a = 100. 



/ 


b 




E\3 31 & (E\3 41) 


hi 7i j 

M 


hi c* r 


n 


i 
i 


op 7n^Q7 


9fi 7n^Qfifi 
zu./uoyuu l^ioj 


op 7n^Qfifi 
zu. i uoyuu 


op 7f\p 
ZU. 1 uu 




1 n 


i n fi 

1U.U 


1 l.OO l ( UJ 


1 1 ^791 Q7 
ii .o 1 z ly i 


11 fi1 1 9 
1 1 .U 1 1Z 




1 nn 

1UU 


o.yozo 


Q Q7Qn (\ 7\ 

o.y i yu \i- 1 ) 


o.yoouyo 






1 nnn 

1UUU 


o.uyyoio 


^ nQQ7QQfi ('sni zL^ 
o.uyy/yyu ^oui^j 


o.uyyoi i 


o.uyyo 


i 
i 


i 
i 


/I9 oi7z>7 


/I9 O^^fil 9 fSS^ 
^tZ.ZOlOUlz ^OOj 


AO 9 f ?7^fin 


/19 OIP 
^tz.ZOU 






1 4 ^7n 


l^i.OUZZ I U^ 1 


1 A ^fiS&l 1 

i'i.OOOOll 


1/1 ^fi^s 




1 nn 

1UU 


^ QQIAIAP 


^ qq'?/! ^nn fen 
o.yyo^touu ^yi ) 


o.yyo^toy 


o.yyou 




1 nnn 

1UUU 


^ nQQQ^^^Zl 


o.uyyyoo^u j 


o.uyyyoo 


^ 1 nnn 

o. 1UUU 


2 


1 


55.97780 


55.9778047 (70) 


55.977804 


55.976 




10 


16.61084 


16.61081 (54) 


16.610869 


16.6110 




100 


7.99602470 


7.99602475 (91) 


7.996025 


7.9960 




1000 


7.099960046 


7.0999600264 (37) 


7.099960 


7.1000 


3 


1 


67.9608094 


67.96080 (81) 


64.960806 


67.9600 




10 


18.719989 


18.720001 (21) 


18.719999 


18.7202 




100 


9.997153662 


9.997153641 (55) 


9.997154 


9.9972 




1000 


9.0999714406 


9.099971440 (39) 


9.099971 


9.1000 


4 


1 


78.2383822 


78.23838049 (45) 


78.238380 


78.2380 




10 


20.7814134 


20.7814161 (67) 


20.781416 


20.7820 




100 


11.9977838 


11.9977838311 (34) 


11.997784 


11.9978 




1000 


11.099977784 


11.099977784 (84) 


11.099978 


11.10000 
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Table 5: Same as table 3 for a = 1000. 



/ 


b 




E\3 31 & (M3 41) 


hi 7i j 

^ M 


hi C T 


n 


i 
i 


Q1 9^fifi 


Q1 9^fifi1 1 1 9 f^9"\ 


Q1 9^fifi1 1 


qi 9^fin 




1 n 


fizL 8^ 




p.A 89^n8^ 
Ul.oZOUOO 


RzL 7zLnn 




1 nn 

1UU 


19 81 R9 


1 9 784 ('en^i 

1Z. 1 O^: ^OU J 


1 9 89^4^ 


1 9 8^fifi 




1 nnn 

1UUU 


Q QQ81 Q 

o.yyoio 


o.yyouu ^ui j 


Q QQ81 n7 
o.yyoiu i 


Q QQ84 


i 
i 


i 
i 


1 /IQ fi^fi'il Q/l 


1/1Q fi^fi'ilQ/IQ ^81 "\ 

i^ty.ooooiy^ty yoi j 


1 /IQ fi^fi^l Q 


1 /IQ fic^n 






8Q 1 9fi 
oy . izu 


8Q 1 9^/LQ (17\ 
oy.izoiy i i 


8Q 1 9 f ?4 c i9 
oy. izo^ioz 


8Q nzinn 
oy .v^w 




1 nn 

1UU 


1 A Q'?'?798 


1 /I Q'i'iQS ('/infi^ 

i^t.yooyo y^vu) 


1 /I 0^^77/1 


i /i Q^n 




1 nnn 

1UUU 


^ QQQ'?'?^'? 

o.yyyoooo 


x> qqq^zl^s (^^^ 
o.yyyoo^oo ^J-J-J 


o.yyyooo 


K QQQ/I 

o.yyy^ 


2 


1 


206.1068055 


206.1068055 (59) 


206.106805 


206.10000 




10 


100.60 


100.7043 (64) 


100.703996 


100.7260 




100 


16.9601067 


16.9601077 (91) 


16.960106 


16.9604 




1000 


7.99960031 


7.999600250 (23) 


7.999600 


7.99960 


3 


1 


260.6091863 


260.6091863 (71) 


260.609186 


260.6200 




10 


105.5097 


105.5054 (70) 


105.507579 


105.4960 




100 


18.9714849 


18.97148472 (84) 


18.971485 


18.9716 




1000 


9.99971440 


9.999714394 (86) 


9.999714 


9.9998 


4 


1 


313.16466655 


313.164666 (72) 


313.164667 


313.1600 




10 


108.52780 


108.527802 (18) 


108.527834 


108.5280 




100 


20.9778139 


20.977813876 (97) 


20.977814 


20.9780 




1000 


11.999777840 


11.999777839 (44) 


11.999778 


11.9998 
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Table 6: k = eigenvalues for V{q) = -l/2q + ([i/2)ln(q 2 + q). 



n 


Ep 


E[4,4] 


STM 


DNI 


0.0001 


0.2497779 


0.2497779 


0.24975 


0.24978 


0.0005 


0.248890 


0.248890 


0.24875 


0.24889 


0.001 


0.247782 


0.247782 


0.24752 


0.24778 


0.005 


0.238976 


0.238973 


0.23765 


0.23897 


0.01 


0.228105 


0.228098 


0.22545 


0.28810 


0.05 


0.145676 


0.145681 


0.13227 


0.14568 


0.1 


0.051363 


0.051499 


0.02456 


0.05153 


0.5 


0.521641 


0.520529 


0.65413 


0.52033 
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Table 7: CUTT OFF POTENTIAL. 



b 




Is 


2p 


3d 


4/ 


0.3 


Ep 


0.2938592647 


0.1061734109 


0.0516188327 


0.0299960986 




£[4,4] 


0.2935835059 


0.1061736042 


0.0516188327 


0.0299960986 




E 

J -'ex 


0.29354528 


0.10617351 


0.05161883 


0.02999610 


0.5 




0.2446947653 


0.0976555780 


0.0515016542 


0.0292431052 






0.2445430901 


0.0976556895 


0.0515016549 


0.0292431052 






0.24453144 


0.09765562 


0.04943696 


0.02924311 


5 




0.0706714094 


0.0434584151 


0.0287056893 


0.0200000000 






0.0706701591 


0.0434584053 


0.0287056891 


0.0200000000 






0.07067028 


0.04345840 


0.02870569 


0.02000000 


100 




0.0067420934 


0.0056337380 


0.0048124934 


0.0041687844 






0.0067420738 


0.0056337379 


0.0048124934 


0.0041687844 






0.00674208 


0.00563374 


0.00481249 


0.00416878 



27 



Table 8: NPO for k = 1 and a = 10. 



b 


I 


EpsLET 


E 44 


El/N 


EcFM 


1 





13.3557146 


13.388417 


13.496 


13.388323 


10 




7.9029106 


7.9654956 


7.9180 


7.903154 


100 




7.0984744 


7.0979219 


7.0990 


7.098449 


1000 




7.0099838 


7.0099799 


7.0100 


7.009982 


1 


1 


16.0066200 


16.0156839 


15.963 


16.016128 


10 




9.9447832 


9.9452031 


9.9504 


9.944898 


100 




9.0993550 


9.0993626 


9.0994 


9.099352 


1000 




9.0099934 


9.0099934 


9.0100 


9.009993 


1 


2 


18.5442860 


18.5432358 


18.498 


18.543473 


10 




11.9633289 


11.9633409 


11.965 


11.963343 


100 




11.0996043 


11.0996041 


11.010 


11.099604 


1000 




11.0099960 


11.0099960 


11.010 


11.009996 
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